Conductance fluctuations in disordered superconductors with broken time-reversal 

symmetry near two dimensions 
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We extend the analysis of the conductance fluctuations in disordered metals by Altshuler, 
Kravtsov, and Lerner (AKL) to disordered superconductors with broken time-reversal symmetry 
in d — (2 + e) dimensions (symmetry classes C and D of Altland and Zirnbauer). Using a per- 
turbative renormalization group analysis of the corresponding non-linear sigma model (NLcrM) we 
compute the anomalous scaling dimensions of the dominant scalar operators with 2s gradients to 
one-loop order. We show that, in analogy with the result of AKL for ordinary, metallic systems 
(Wigner-Dyson classes), an infinite number of high-gradient operators would become relevant (in 
the renormalization group sense) near two dimensions if contributions beyond one-loop order are 
ignored. We explore the possibility to compare, in symmetry class D, the e = (2 — d) expansion 
in d < 2 with exact results in one dimension. The method we use to perform the one-loop renor- 
malization analysis is valid for general symmetric spaces of KaMer type, and suggests that this is a 
generic property of the perturbative treatment of NLctMs defined on Riemannian symmetric target 
spaces. 



I. INTRODUCTION 

A. Conductance fluctuations and the non-linear 
sigma model 

The non-linear sigma model (NLcrM) provides, as first 
proposed by F. Wegnerji an efficient framework within 
which to formulate (Anderson-type) localization tran- 
sitions of noninteracting quantum mechanical particles 
subject to (static) random potentials. This provides some 
of the simplest models of electrons in disordered solids. In 
a seminal piece of work, Altshuler, Kravtsov and Lerner 
(AKL)-i 3 . showed in the 1980's how to incorporate the 
important phenomenon of conductance fluctuations^^ 
into this framework. Within a perturbative renormaliza- 
tion group (RG) analysis of the weakly coupled NLcrM in 
d = 2 (or d — 2 + e) spatial dimensions, i.e., in the regime 
of large (dimcnsionlcss) conductance g, these authors 
stressed the importance of terms in the low-energy effec- 
tive action of the NLcrM which possess an arbitrary large 
number (= 2s) of spatial derivatives but do not violate 
any symmetries. These terms are called "high-gradient 
operators" and are generalizations of the kinetic term in 
the NLcrM which has two derivatives (2s = 2). High- 
gradient operators are highly irrelevant in the RG sense 
by "power counting" (i.e., without incorporation of fluc- 
tuation corrections). Specifically, AKL and Yudson stud- 
ied the anomalous (RG) dimensions of these operators to 
one-loop order in the fluctuation corrections in d = 2 and 
in the e = (d — 2)-cxpansion for d > 2£&Z& It was found 
that (i): these operators would become strongly relevant 
(in the RG sense) for a large enough number 2s of deriva- 
tives, in spite of their strong irrelevance by power count- 
ing, if contributions beyond one-loop order are ignored, 
and that (ii): they would (consequently) dominate the 



cumulants of the conductance of sufficiently large order 
2s, in a sample of (large) linear size L. Based on these 
observations AKL were led to argue that for large mean 
(dimcnsionlcss) conductance g, the main weight of the 
probability distribution of the conductance arises from a 
narrow peak at g = g, while there exist anomalously long 
tails (close to being log- normal). The latter arise from 
the high moments of the conductance, which are dom- 
inated by the above-mentioned high-gradient operators 
appearing in the low-energy effective action. Upon ap- 
proaching the regime of stronger disorder (smaller mean 
conductance g), and in particular upon approaching a 
metal-insulator transition, the probability distribution 
becomes increasingly broad. AKL find similar behavior 
for the probability distributions of other obscrvablcs, in- 
cluding, e.g., the local density of states (LDOS), current 
densities, and long-time tails of relaxation currents. Very 
recently^S the basic observations of AKL for the LDOS 
were confirmed in a chiral symmetry class (nomenclature 
of Rcfs. [ll] and|T3), and were analyzed more completely 
in this case. 

Ten years ago, in another seminal piece of work, 
Zirnbauerjii and Altland and Zirnbauefti 2 . (AZ) pro- 
vided a general classification scheme for the behavior 
of non-interacting quantum mechanical particles sub- 
ject to random potentials, on length scales much larger 
than the mean free path in terms of a total of ten 
symmetry classes*^ In every realization of disorder, the 
Hamiltonian has the same symmetries as a correspond- 
ing ensemble of random matrices. This encompasses 
the three previously known so-called "Wigner-Dyson 
symmetry classes" relevant for the physics of Ander- 
son localization of electrons in disordered solids (corre- 
sponding to "orthogonal, unitary, and symplectic" ran- 
dom matrix ensembles)^ Altland and Zirnbauer's exten- 



2 



sion includes four novel symmetry classes describing the 
(Anderson-like) localization physics of non-interacting 
Bogoliubov-De Gennes (BdG) quasiparticlcs within a 
mean-field treatment of pairing in disordered supercon- 
ductors. (Four symmetry classes arise since SU(2) spin- 
rotational, or time-reversal invariance may be present or 
absent.) As parameters are varied, the wave functions 
of the (BdG) quasiparticlcs may be extended or local- 
ized due to disorder, and (Anderson localization type) 
transitions between such phases may occur. Since BdG 
quasiparticles do not carry a conserved electric charge, 
the difference between phases of localized and extended 
quasiparticle wave functions manifests itself not in the 
electrical conductivity, but instead in the thermal con- 
ductivity k, and if spin is conserved, also in the spin 
conductivity er spin J£ In complete analogy with the fluc- 
tuations of the electrical conductance in a finite sam- 
ple of disordered metals (described by the Wigner-Dyson 
symmetry classes) , the system of BdG quasiparticles also 
exhibits fluctuations of the appropriate conductance in 
a finite-size sample. These arc the fluctuations of the 
thermal conductance (divided by temperature T), or of 
the spin conductance. Henceforth, for notational con- 
venience, we will refer to these quantities simply as "the 
conductance" , and denote the corresponding "dimension- 
less conductance" again by the symbol g. 

Just as for disordered metals (i.e., for the Wigner- 
Dyson symmetry classes), the NLcrM also provides an 
efficient framework to formulate (Andcrson-typc) local- 
ization physics of BdG quasiparticles in superconductors. 
In this article we focus on systems where time-reversal 
symmetry is broken in every realization of disorder: the 
corresponding AZ symmetry classes are conventionally 
called "class D" if SU(2) spin-rotational invariance is 
broken and "class C" if it is conserved, in every disor- 
der realization^ In parallel, and for comparison, we will 
also discuss here localization physics of ordinary elec- 
trons in the absence of time-reversal symmetry, which 
is described by the Wigner-Dyson symmetry class con- 
ventionally called "unitary class" or "class A" . (See also 
Table fl] ) When disorder averaging is implemented using 
fermionic replicas, the target manifolds of the resulting 
NLerMs are the symmetric spaces 0(2iV)/U(iV) for sym- 
metry class Djii and Sp(7V)/U(7V) for symmetry class 
CJ^ We remind the reader that for the familiar unitary 
symmetry class A, the corresponding target spaced is 
U(2A)/[U(A0 x U(A)]. (In all three cases the number 
N of replicas is taken to zero at the end of calculations.) 

The aim of this paper is to investigate perturbativcly, 
in the spirit of the original work of AKL for ordinary 
metals, the conductance fluctuations for the novel AZ 
symmetry classes D and C, i.e., for BdG quasiparticles in 
superconductors lacking time-reversal symmetry, in the 
absence (class D) or presence (class C) of spin-rotational 
symmetry. Our main results are briefly summarized in 
Sec. II B I below. We also point out that it might be possi- 
ble in the case of symmetry class D to make contact with 
exact nonperturbative results obtained in d = 1 dimcn- 



RMT T SU(2) 


NLcrM target space 


■d 


TV 


D no no 


O(27V)/U(A0 


+1 


N 


C no yes 


Sp(JV)/U(JV) 


-1 


N 


A no 


U(p + g )/[U(p)xU( g )j 





p + q 



TABLE I: Fermionic replica target spaces of the NLcrM de- 
scribing weakly disordered superconductors and metals when 
time-reversal symmetry (T) is broken. Symmetry classes D, 
C, and A refer to the symmetric spaces defined by the RMT 
limits of the NLcrM. The number of replicas M is N for the 
two BdG symmetry classes whereas it is p = q — N for the 
unitary symmetry class. The index $ = 0,±1 distinguishes 
the behavior of the beta function for the NLcrM in symmetry 
classes D, C, and A. 

sions from the so-called Dorokhov-Mcllo-Pereyra-Kumar 
(Fokkcr-Planck) analysis — of the probability distribution 
of the conductance. 

We end this section by recalling, as an aside, that the 
NLerMs in all three symmetry classes, D, C, and A, share 
certain essential features associated with time-reversal 
symmetry breaking. First, in d = 2 dimensions, all three 
target manifolds allow for a topological term (a term that 
does not modify the equations of motion), the so-called 
"theta- or Pruisken term". Its presence is responsible 
for the appearance of thermal, spin, or charge Hall insu- 
lating phases, and corresponding plateau transitions in 
d = 2 dimensionsi 18 i 21 i 22 ' 23 i 24 Second, the target mani- 
fold for symmetry class D stands out in that it is not 
simply connected. 11 ' 25 ' 26 Consequently, there are certain 
discrete Z 2 domain wall excitations that may give rise to 
a rich phase diagram for certain types of disorder. 23 ' 27 ' 28 
In this paper we shall ignore the physics associated with 
such Z 2 domain walls, which requires an extension of 
the NLcrM description in symmetry class D. These fea- 
tures are not present in certain realizations of this sym- 
metry class, on which we focus in this paper. Finally, all 
three target spaces have in common that they arc Kahler 
manifolds from a geometrical point of view. This means 
that each target space defines a complex manifold en- 
dowed with a Hcrmitian metric derived from a Kahler 
potentia h 29 ' 30 This property will allow us to perform the 
perturbative RG analysis in a unified way for symmetry 
classes D, C, and A. (See Sec. IIP I for a brief summary.) 

B. Summary of main results 

We begin with the summary of our main results by 
quoting the standard one-loop resul t 32 ' 36 for the (RG) 
beta function^ 3 . 

0{t) = -et + {AT - d)t 2 + 0{t 3 ) (1.1) 

in d = 2 + e spatial dimensions, obeyed by the dimen- 
sionlcss coupling constant t of the relevant NLcrM [de- 
fined in (|2.1b[) below] whose inverse is proportional to 
the mean dimcnsionlcss conductance g, once the replica 
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Class C, d > 2 Class C, d < 2 

t t 

gm — « — . m - — - — *- g— — > — > — > — <■ — - — 

t* 

Class D, d > 2 Class D, d < 2 

t t 
G— « « — « &* > m ' — «- 

t* 

FIG. 1: Infrared (IR) renormalization group (RG) flows for 
symmetry classes C and D after taking the replica limit de- 
pending on whether dimensionality of space d is larger or 
smaller than 2. The open circle at t = is the metallic fixed 
point. The filled circle at t* is the metal-insulator critical 
point in symmetry class C whereas it describes a stable fixed 
point for symmetry class D. Up to one loop, t* — |e| where 
e := d — 2 is positive for symmetry class C while it is negative 
for symmetry class D. 



d 




FIG. 2: Possible evolution of the stable fixed point for sym- 
metry class Dind=2 — |e| < 2 dimensions between d = 1 
and d = 2 (after taking the replica limit). In this sketch it 
is assumed that the zero t* = |e| + C(e 2 ) of the beta func- 
tion, found to one-loop order in the e expansion in d — 2 — |e| 
dimensions, defines a curve of stable critical points which in- 
terpolates between the metallic phase in d = 2 dimension (at 
zero coupling t = 0), and the behavior of the corresponding 
symmetry class D theory in d = 1 dimension, which can be 
solved exactly. 



limit TV — > has been taken. To treat the symmetry 
classes D, C, and A on equal footing, we have introduced 
the replica index M that takes the value N for symme- 
try classes D and C and p + q for symmetry class A. 
The index $ is also needed to distinguish the one-loop 
beta function for the NLerM corresponding to symmetry 
classes D (d = +1), symmetry classes C (d = —1), and 
symmetry classes A ($ = 0). The correspondence be- 
tween the pair ($,AT) and the symmetry classes D, C, 
and A is summarized in Table [J The beta function in 
Eq. (|1.1[) possesses a trivial fixed point, 



t = 0, 



(1.2a) 



which is infrared (IR) stable for d = 2 + e > 2 (i.e., 
e > 0), while it is IR unstable when d = 2 + e < 2 
(i.e., e < 0). Upon taking the replica limit Af — ► 
this describes a diffusive metallic phase. Whenever < 
e/(Af — i?) < oo, there is a nontrivial zero t* of the beta 
function (fl~T|^ ^^^ at 



0{e 2 



;i.2b) 



For symmetry classes C and A this describes in the 
replica limit and inrf = 2 + e>2 dimensions, a metal- 
insulator transition separating a metallic from an insu- 
lating phase. For symmetry class D, on the other hand, 
this describes (in the replica limit) a stable fixed point 
("critical phase") in d = 2 — |e| < 2 dimensions. The IR 
flows for symmetry classes C and D are summarized in 
Fig. [1] A possible evolution of the stable fixed point in 
symmetry class D for dimensionality d = (2 — |e|) between 
one and two dimensions, to be discussed in the following 
subsection, is sketched in Fig. 

In this paper we find the dominant anomalous scaling 
dimension among all high-gradient operators con- 



taining 2s gradients, to one-loop order. These are 



2s-2(s 2 -s)t + 0(t 2 ) 

forU(p + g)/[U(p) xU(?)l, (1.3a) 



2s-2(s 2 -s)t + 0{t 2 ) 
for Sp(iV)/U(iV), 
(s) = 2s- (s 2 -s)t + 0(t 2 ) 

for SO(2A)/U(A), 



(1.3b) 
(1.3c) 



when evaluated in the vicinity of the zero-coupling fixed 
point (|1.2a[) where t is small. At the corresponding fixed 
points in d = 2 + e dimensions, Eq. (jl.2bp . these anoma- 
lous scaling dimensions become, consequently [for generic 
numbers of replicas, (p, q) or N] 



2s - 2(s 2 - s)- 



0(e 2 



p + q 

forU(p + g)/[U(p) xU( 9 )l, (1.4a) 



2s - 2(s 2 - s) 



N + l 
for Sp(A)/U(A), 



+ 0(e 2 ) 



2s - (s 2 - s) 



0(e 2 



N - 1 
for SO(2A)/U(A), 



(1.4b) 
(1.4c) 



Equations (|1.3a|) and (|1.4a|) were derived a long time ago 
by Lerner and Wegner . 38 ' 39 Equations (|1.3p and (|1.4p are, 
otherwise, new and are the main results of this paper. 

Taking the replica limit of Eq. (|1.4p yields the dom- 
inant anomalous scaling dimensions at the nontrivial 
fixed points of the localization problems in d = 2 + e 
dimensions^ In particular, in symmetry class C one ob- 
tains, in the replica limit, the one- loop anomalous scaling 
dimensions 



) = 2s-2(s 2 -s)e + 0(e 2 ) 



(1.5a) 
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of dominant 2s-gradicnt operators at the metal-insulator 
transition in d = 2 + e > 2 dimensions (e > 0). On 
the other hand, in symmetry class D in ef = 2 — |e| < 2 
dimensions, the dominant 2s-gradient operators have the 
anomalous scaling dimensions 

x& = 2s-(s 2 -s)|e| + 0(e 2 ) (1.5b) 

at the stable, non-trivial fixed point (in the replica limit). 

C. Discussion 

When taken at face value, the one-loop expressions in 
Eqs. p.3p . (|1.4p . and p.5p would appear to imply that 
the anomalous scaling dimension x^ always turns neg- 
ative for sufficiently large s4£ If this behavior was to 
persist when all higher orders in e are included, an in- 
finite number of high-gradient operators would be rele- 
vant at the nontrivial fixed points of these NLcrMs, i.e., 
infinitely many coupling constants would have to be fine- 
tuned to reach these nontrivial fixed points. The per- 
turbative treatment of these NLcrMs would then appear 
to contradict the belief that the underlying microscopic 
models undergo continuous phase transitions whose ap- 
proach obeys single-parameter scaling laws. Of course, 
higher order contributions in e to these anomalous scal- 
ing dimensions, not computed in the present or earlier 
work, may well change this property and could render the 
actual anomalous scaling dimensions x^*' larger than the 
spatial dimension d — 2 + e, thereby making them irrele- 
vant^-—————— [Indeed, similar behavior is known 

to occur in other random disordered systems (see e.g. 
Refs. E2 and[I3). - See also the discussion in Ref. HI ] 

Nonperturbative properties of the metal-insulator 
transitions in d = 3 dimensions have been obtained 
numerically for the most part . 48 ' 49 i 50 ' 51 They support 
the one-parameter scaling hypothesis for metal-insulator 
transitions in d — 3 dimensions, whose critical behav- 
iors are presumed to be described by the non-trivial fixed 
point of the replica limit of the NLerM in the correspond- 
ing symmetry classes. Whether the apparent breakdown 
of one-parameter scaling due to the relevance of high- 
gradient operators at one loop is an artifact of the 2 + e 
expansion or has a deeper meaning remains an outstand- 
ing problem for the description of Anderson localization 
using the NLerM approach. 

Symmetry classes D and Dili, on the other hand, may 
perhaps offer a somewhat different opportunity for the 
following reason. On the one hand, there exists a critical 
phase described by a stable fixed point in d = 2 — |e| < 
2 dimensions whose properties are accessible perturba- 
tively within the e expansion. These include the prop- 
erties of the conductance fluctuations and of the high- 
gradient operators discussed above. On the other hand, 
a critical phase of these models is also known to exist in 
d = 1 dimensions. Therefore, one might expect that the 
results obtained from the e expansion in d = 2 — |e| < 2 
dimensions would eventually, when |e| = 1, turn into the 



results obtained from the exact solutions of the models in 
d = 1 dimensions, as is schematically depicted in Fig. [5] 
If this was the case, then one could "test" the proper- 
ties of the high-gradient operators and the nature of the 
conductance fluctuations obtained from them following 
the work of AKL, against exact results in d = 1 dimen- 
sions. Of course, as with all e expansions, one would 
not be able to predict exactly the actual numerical val- 
ues of critical exponents directly at |e| = 1. However, 
based on experience with the e expansion in, say, ordi- 
nary Landau-Ginzburg theory, one might expect that the 
structure of the entire probability distribution of the con- 
ductance would evolve, in some sense, "smoothly" with e, 
possibly all the way down to d = 1 dimensions (|e| = 1). 
Let us now pursue this possibility. 

Using perturbation theory of the NLerM, AKL demon- 
strate that in general (i.e., for any NLerM) there are two 
contributions to the sth cumulant of the (dimensionlcss) 
conductance^ (Below, we write the corresponding ex- 
pressions appearing directly at the stable fixed point in 
d = 2 — |e| dimensions.) The first, "standard" contribu- 
tion reads^ 2 - 

«.9 s » s td cx {gf- s (1.6) 

where g = (g) is the mean conductance. However, when 
a high-gradient operator with 2s derivatives is added to 
the action of the NLerM with coupling constant z s , then 
there is an additional contribution 

«S*»add « Zs (jj (1-7) 

to this cumulant. Here, L is the linear system size and I 
the characteristic microscopic length scale (the "mean 
free path") for the disorder. We see that it is when 
d — x^ s > > that the additional contribution in Eq. (|1.7p 
dominates over the standard contribution in Eq. (|1.6[) . 
This hence happens precisely when the high-gradient op- 
erator with 2s derivatives becomes a relevant operator in 
the action of the NLerM. 

Now, AKL argue that, if one was to take the one-loop 
expressions for the anomalous scaling dimensions of the 
high-gradient operators at face value, which for symme- 
try class D would amount to the result in Eq. (|1.5b() 
derived by us in this paper, then for cumulants of suf- 
ficiently higher order s > jfr[l + O(e)] the additional 
contribution would dominate over the "standard" one, 
yielding 

((9 S ))-A s ^y X , s> + (1-8) 

as L/£ — > oo, where A s are certain amplitudes. Following 
AKL, the lower cumulants of order s < ||y[l + 0(e)], on 

the other hand, would be determined by Eq. (jl.6p . In 
summary, the reasoning of AKL implies, when specialized 
to the stable fixed point in symmetry class D in d = 
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2 — |e| < 2 spatial dimensions, the following form of the 
conductance cumulants 

f {g) 2 - s , if s < so, 

((9 s )) ~ { (L9a) 
{ A 8 (L/£) d ~ x , if s>s , 

where 

a0 «A[l + O( c )] (1.9b) 

and 

d - z« = |e|s 2 - (2 + |e|)s + (2 - |e|) + 0(e 2 ). (1.9c) 

Let us now discuss the known exact results for the 
d = 1 dimensional case. Specifically, the leading depen- 
dence of all cumulants of the conductance of a very long 
but thick quasi one-dimensional wire of length L, be- 
lieved to be described by a NLcrM in d = 1 dimensions, 
has been extracted exactly from the DMPK approach for 
symmetry classes D and DII L 53 i 54 i 55 For both symmetry 
classes, these cumulants decay algebraically with system 
size, 

((<?*)) oc (L/i)- 1 / 2 , 8 = 1,2,..., (1.10) 

with the same exponent for all cumulants ("£" is again 
the mean free path). Moreover, the logarithm of the 
conductance exhibits the behavior 

((In ff)) cx -(L/£) 1/2 , vaxhig~L/l (1.11) 

Thus, the typical conductance exhibits stretched expo- 
nential behavior, in contrast with the mean conduc- 
tance. This indicates a broad conductance distribution. 
[For symmetry class Dili the result (| 1 . 10[) for the mean 
conductance has been reproduced by Lamacraft, Si- 
mons, and Zirnbauer^ using the one-dimensional NLcrM 
(which is a principal chiral model for symmetry class 
Dili) and instanton methods.] 

If one assumes that the exact results in d = 1 describe 
the end point (at e = 1) of the line of critical points in 
2 — |e| dimensions (see Figj2j) , there are then two possibil- 
ities. First, the one-loop relevance of high-gradient oper- 
ators is an artifact of the e-expansion and these operators 
are in fact irrelevant once all orders in the e-expansion 
are taken into account (see Refs. [jiliiil ). Second, this is 
not the case and high-gradient operators are truly rele- 
vant. In this case, however, deducing the scaling of the 
cumulants based on the one-loop anomalous scaling di- 
mensions does not appear to be sufficient to determine 
their scaling behavior. Indeed, in this case the running 
coupling constants z s of all the high-gradient operators 
in the effective action of the NLcrM grow rapidly, more 
so the higher the number 2s of gradients. Under these 
conditions, then, there is no rationale to neglect the non- 
linear coupling between the running coupling constants 
z s of the high-gradient operators as was done to reach 



the conclusion (|1.9[) . when extracting the scaling of the 
cumulants of the conductance. In effect, one needs to 
implement a functional RG analysis which can treat an 
infinite set of running (and increasing) coupling constants 
on equal footing. We have performed such an analysis (on 
which wc will report elsewhere) for symmetry class A, but 
have so-far been unable to solve the resulting RG equa- 
tions. Since in the present case the coupling constants 
z s of the high-gradient operators all become large, it is 
clear (as already mentioned) that they will not flow inde- 
pendently of each other in the IR limit, but will strongly 
mix with each other under the RG flow. In this situa- 
tion, their flow does not need to result in the breakdown 
of single parameter scaling. Indeed, a precisely analo- 
gous situation is encountered when computing the local 
DOS of the two-dimensional random hopping problem^ 
in which case one-parameter scaling for the local DOS is 
recovered using a functional RG analysis. 



D. Summary of the method of calculation 

The RG scheme that we adopt in this paper is the co- 
variant background field method^ In Refs. HI45I this 
method was used to compute, up to two loop order, 
the RG flows of the high-gradient operators for the 
0(N)/0(N - 1) NLcrM, relevant to statistical models 
such as the classical Heisenberg magnet. The covariant 
background field method consists of two key steps. First, 
the fields of the NLcrM are separated into short wave- 
length ("fast") and long wavelength ("slow") degrees of 
freedom where the slow degrees of freedom are defined as 
the solutions to the classical equations of motion. Sec- 
ond, the fast fields are integrated out. When using this 
method, no redundant operator— is generated in the sec- 
ond step of the covariant background field method, a great 
advantage when one needs to consider mixing of many 
scaling operators by the RG transformation (as we do 
here). Another advantage is that explicit invariance un- 
der a coordinate reparametrization of the group manifold 
is maintained, allowing for a geometrical formulation at 
all stages of the calculation. 

In this paper, we shall also extend the covariant 
background field method to the cases when the tar- 
get spaces of the NLcrM are Ricmannian or Kahler 
supermanifoldsi^ We derive a master formula, Eq. (|4.3[) . 
in the Kahler case for the RG flow of high-gradient opera- 
tors, all of which possess a vanishing covariant derivative. 
[See Eq. (|C19|) for the super Riemannian^ and Eq. (|C40|) 
for the super Kahler— cases, respectively. To obtain the 
corresponding formula for the non-supersymmetric cases 
(for fermionic replicated NLcrM, i.e., for compact NLcrM) 
we set all factors (— l) p to +1, where p = 0, 1 is here the 
grade in the supersymmetric expression]^ This master 
formula is expressed solely in terms of the geometrical 
data of the target manifold and should apply to all 10 
symmetry classes of AZ ; 11 ! 12 but is not limited to them. 
The generality of the master formula suggests that high- 
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gradient operators are generically relevant to 1-loop or- 
der at the nontrivial fixed point of a NLcrM defined on a 
Ricmannian target space. 

E. Geometric interpretation of high-gradient 
operators 

All high-gradient operators that enter the computation 
of the cumulants for the conductance can be expressed 
in terms of tensor fields on the target manifold. We will 
show that high-gradient operators together with the met- 
ric and curvature tensors on the target manifolds of the 
NLcrM form a family of tensors with vanishing covariant 
derivative. 



F. Outline of the paper 

This paper is organized as follows. Section |TT] intro- 
duces the NLcrM on Kahler manifolds for disordered su- 
perconductors with broken time-reversal symmetry as 
well as for a disordered metal in the standard unitary 
symmetry class. Section IIIII describes how to use the 
covariant background field method to deduce pcrturba- 
tive RG flows. The dominant anomalous scaling dimen- 
sions for high-gradient operators are computed in Scc. lIVI 
to one-loop order. We relegate to appendices some key 
intermediary steps in the derivation of Eq. (|1.3p . The 
master formula, Eq. (|4.3|) . which expresses the RG trans- 
formation law obeyed by high-gradient operators solely 
in terms of the geometrical data of the target manifold is 
derived in appendix [C] for any Riemannian or Kahlcrian 
supermanifold on which the high-gradient operators have 
vanishing covariant derivatives. 

II. NLcrM FOR WEAKLY DISORDERED 
SUPERCONDUCTORS 

A. NLcrM on Kahler manifolds 

We start from the NLcrM description for the symme- 
try classes D, C and A of Anderson localization in terms 
of fermionic replicas. Alternatively, one can use either 
bosonic replica or supersymmetric methods to perform 
the disorder average. These three options are equivalent 
as far as perturbative RG computations are concerned. 
The field entering the NLcrM is a map from the base 
space, a d-dimensional Euclidean space M. d , to the target 
manifold 9Jt defined by the pattern of symmetry break- 
ing. The target manifold has a well-defined natural met- 
ric g. A classification of all the target manifolds of NLcrM 
relevant to the problem of Anderson localization applied 
to non-interacting quasi-particlcs were identified in the 
seminal work of Zirnbauer, and Altland and Zirnbauer.— 
For the cases of interest in this paper, i.e., disordered su- 
perconductors with broken time-reversal symmetry, the 



target spaces of the relevant NLcrM turn out to be Kahler 
manifolds which are the Hermitian symmetric spaces 
SO(2A0/U(A0 for symmetry class D and Sp(N)/\J(N) 
for symmetry class C^- The target space of the NLcrM 
that describes a disordered metal when time-reversal 
symmetry is broken is the Hermitian symmetric space 
U{2N)/[\J(N) x U(A r )], which is isomorphic to the com- 
plex Grassmannian G 2N N (C)£Q 

Within the mesoscopic community, the usual represen- 
tation of the NLcrM in symmetry classes D, C, and A 
is in terms of matrices belonging to the relevant sym- 
metric spaces, "the Q-matrix representation" after the 
seminal work by Efetov et al. in Ref. 0. This repre- 
sentation emphasizes the pattern of symmetry breaking 
but quickly becomes cumbersome when performing an 
RG analysis of the flow of composite operators. We opt 
instead to emphasize more directly the geometrical prop- 
erties of the target spaces, namely the fact that they are 
examples of Kahler manifolds for symmetry classes D, 
C, and A. We refer the reader to appendix [A] for a re- 
lation between the representation of the NLcrM that we 
shall introduce below and the more common Q-matrix 
representation thereof. 

We now present the definition of a NLcrM on a Kahler 
manifold starting from a NLcrM on a 2M-dimcnsional 
Ricmannian manifol d 30 ' 31 (9Jl,g), 

Z := J V[cf\e- S[ *\ (2.1a) 

m := 9A B {4>)d^ B d^ A . (2.1b) 

Our conventions are here the following. We reserve 
the greek alphabet to label the d coordinates r where 
/i = 1, . . . , d of r 6 M. d , and use the short-hand notation 
J r = J| r | >0 d r/a 2 for the integration over R d with 
the short-distance cut-off a. The dimensionless coupling 
constant of the NLcrM is here denoted by t. The dimen- 
sionless conductance shares with (<\ d ~ 2 t)~ l the same en- 
gineering dimension. To each point r from the base space 
we have assigned a point p on the manifold 971 with the 
coordinates (<f> A ) G R 2M . In the definition of the NLcrM 
on a Ricmannian manifold, Eq. (|2.ip . the 2M x 2M com- 
ponents of the metric tensor field g of rank (0, 2) are rep- 
resented by (g AB )- We reserve the capital latin alphabet 
to label the 2M coordinates ^=l.-.2M f p e 971. 

The components of the metric tensor field g AB on the 
Riemannian manifold (971, g) are real and symmetric un- 
der the exchange of the indices 

9ab=9ba, A,B = 1,...,2M. (2.1c) 

There is a distinction between upper/lower capital latin 
indices. The components of are denoted by g AB and 
satisfy 

9ac9 CB =aS B , 9 AC 9cb = A 5b, (2-2) 

and so on for A,B,C = 1 . . . , 2M. On the other hand, 
there is no distinction between upper and lower greek in- 
dices and we will always choose them to be lower indices. 
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A Hcrmitian manifold (SEJl, 9) is a 2M-dimcnsional Rlc- 
mannian manifold endowed with a rank (1,1) tensor field 
Z that is defined globally on 971 and is called a complex 
structure. On a Hermitian manifold, we can always make 
a coordinate transformation from (<p A ) to (x a ,y a ) that 
satisfies, for each point p £ 97t, 



i) 



d 



d 



d 

dx a ' 



(2.3) 



where a = 1, . . . , M. It is then natural to introduce com- 
plex coordinates (z* a ,z a ) through 

z* a := x a - iy a , z a :=x a + iy a , a = 1, . . . , M. (2.4) 

A NLctM on a Hermitian manifold (97t,fl) is thus defined 
by the partition function 



Z := 



V\z\, 



-S[z*,z] 



(2.5a) 



51 



**,*] ~ 5a- h (**>*) W°> ( 2 - 5b ) 



where we have assigned to each point r from the base 
space a point p on the manifold 971 with the coordinates 
(z* a ,z a ) e C M . Components of tensor fields on 97t have 
now two types of indices, holomorphic (a, 6, . . .) and an- 
tiholomorphic ones (a*, 6*, . . .). Holomorphic indices can 
be contracted with holomorphic components of the co- 
ordinates z a whereas antiholomorhic ones with z* a , the 
complex conjugate of z a . In the definition of the NLctM 
on a Hcrmitian manifold, Eq. (|2.5[) . the M x M compo- 
nents of the Hermitian metric tensor field g of rank (0, 2) 
are represented by (g a » b ). We reserve the latin alphabet 
to label the M holomorphic and antiholomorphic coor- 
dinates z* a and z a where a*, a = 1, ... , M, respectively, 
of p G 971. The use of a capital latin letter on a Her- 
mitian manifold, say A, as an index denotes either an 
holomorphic (A = a) or an antiholomorphic (A = a*) 
index. 

The components of the metric tensor field g AB on the 
Hermitian manifold (SEJl, q) are symmetric under the ex- 
change of the indices and Hcrmitian 



9ab — 9bai (9ab) — 9a*b* 



(2.5c) 



with A = a or a*, B = b or b*, A* = a* or a, and B* = b* 
or b running from 1 to M. The components of g -1 are 
denoted by g AB and satisfy 



9ac*9 C ' B =a5 B , 9 AC *9c*b= A 5 



B ■ 



(2.6) 



and so on for A, B, C* = 1, . . . , M. 

A Kahler manifold (971, g) is a Hermitian manifold 
whose Hermitian metric can be derived from a Kahler 
potential K(z*, z), 



1ab-( z *, z ) = d a d b ,K{z*,z), 



(2.7) 



with a, b* = 1 . . . , M. Not all Hermitian manifolds are 
Kahlcrian. 



The Kahler manifolds corresponding to the Hermi- 
tian symmetric spaces \J(p + q)/{XJ(p) x U(g)] for sym- 
metry class A, Sp(A)/U(A) for symmetry class C, 
and SO(2A)/U(A) for symmetry class D are specified 
once the corresponding Kahler potentials have been con- 
structed. To this end, we begin by parametrizing the 
complex Grassmannian U(p + g)/[U(p) x V(q)} for sym- 
metry class A by the stereographic coordinates defined 
by the independent set 

W*Aa^Aa, A=l t ...,p, d=l,...,q} (2.8a) 

of p x q complex- valued matricesj 62 i 63 ' 64 With this coor- 
dinate system, the Kahler potential on \J(p + q)/[U(p) x 
U(q)] is given by 

K(<p*,tp) = clndct(l g + (^V) (2.8b) 

where c is an arbitrary constant (the scale of the metric) 
which is chosen to be unity in the following. (Another 
choice for c just amounts to a rescaling of t.) The Kahler 
metric induced by the Kahler potential (|2.8[) is 

S(Aa)(B6)«(¥'* ! ¥>) ; 

where we have introduced q x 
matrices, 



Zab^BA 



(2.9a) 
and p x p Hermitian 



Y 



AB 



ab ' 

•ab 



a, b = 1, , 
A,B = 1,. 



,q, (2.9b) 
,p. (2.9c) 



The cases of Sp(A)/U(A) for symmetry class C and 
SO(2A)/U(A) for symmetry class D are obtained by set- 
ting p = q = N and by imposing the constraint 



ip + -dip 1 







(2.10) 



on the complex Grassmannian XJ(p + q)/[\J(p) x V(q)} 
where $ = — 1 for Sp(A) (symmetry class C) and = 
+1 for SO(2A) (symmetry class D). Due to the con- 
straint (|2.10p . independent entries of {^>\ a , <f 'A a , A, a — 
1,...,N} can be chosen to be those with the indices 
A < a for Sp(A) (symmetry class C) and A < a for 
SO(2A) (symmetry class D), say. In turn, the metric 
for Sp(A)/U(A) (symmetry class C) and the metric for 
SO(2A)/U(A") (symmetry class D) are given by 



a Aa a Bb^(Aa){Bb)^ab^ABy 



where 



°(Aa)(Bb)--- 

^fl rp 



= 1-W2, 

— °(Aa)°(Bb) 



T 



■■■A- 



1?T 



(2.11a) 

(2.11b) 
(2.11c) 

(2.11d) 



and 



Z ab := (l N + <p^) ab , a, b = 1, . . . , N. (2.11e) 

Observe that Z in Eq. (|2.9b|) and Y in Eq. (|23c]) arc not 
independent anymore because Y = Z T when p = q = N 
and the constraint (|2.10p holds. By a slight abuse of 
notation we will sometimes refer to the complex Grass- 
mannian case by i9 = 0. 
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B. High-gradient operators 

The seminal contribution by Altshuler, Kravtsov, and 
Lerner^ was to observe that the perturbative RG flow 
obeyed by the cumulants of the conductance necessarily 
involves an infinite number of composite operators and 
to identify a criterion for the selection of this infinite 
family of composite operators. The selection criterion is 
that the flow must be closed within a family of compos- 
ite operators characterized by the most dominant scaling 
dimensions. Applied to the three cases at hand, we thus 
seek a family of composite operators that transform as 
scalars under reparamctrization of the target manifold 
and contain as many derivatives as possible subject to 
the condition that their anomalous scaling dimension be 
dominant. We are led to the family of high-gradient op- 
erators 



T IJK „Xz*,z)dz I Mdzf--; 



(2.12) 



2.s 



where T IJK (z*,z) is a tensor on the target manifold 
with s holomorphic and antiholomorphic indices; we do 
not yet specify how the greek indices must be contracted. 
Although composite operators containing derivatives of 
the form d n z (n > 2) are generated by the RG flow 
for the cumulants of the conductance they are neither 
expected to contribute to the dominant RG anomalous 
scaling dimensions^ nor needed to close the RG flows of 
the high-gradient operators (|2.12p as we shall sec below. 

As the simplest example of a high-gradient opera- 
tor belonging to the family (|2.12p is the Lagrangian 
g IJ d fl z J d^z I of the NLerM action. This is in fact the only 
scalar under both reparamctrization of the target man- 
ifold and global rotations in the base space with s = 1 
that belongs to the family (|2.12[) . To see this, observe 
that the number of independent tensors with two indices 
on the target manifold is related to that of independent 
bilinear forms in the group (cosct). For simple Lie groups 
(and their analogues for cosets), there is only one such 
bilinear form (the Casimir), and hence gjj is uniquely 
fixed. The scaling flow obeyed by g IJ d )1 z J d ll z I , controls 
the scaling flow obeyed by the mean conductance, i.e., 
the first cumulant of the conductance once the replica 
limit has been taken. 

Life already becomes more complicated if we seek all 
the high-gradient operators with ,s = 2 in the fam- 
ily (|2.12|) . On a generic complex manifold more than 
2 independent tensors of rank (0, 4) are available, 



9lj9KLy Ri,jkl-> ° 



IJMN 



R 



MN 



KLi 



(2.13) 



and so on, where Rjjxl IS the curvature (Ricmami) ten- 
sor on the target manifold. Fortunately, the situation 
is simpler for the Grassmannian since we only need to 
contend with two independent tensors of rank (0,4), 



S 



(i) 



(2) 

(Aa)(Bb)(Cc)(Dd) 



,5' 



(i) 



(Aa)(Bb) , -'(Cc)(Dd) 



(2.14a) 
Z ab Y BA Z cd Y DC . (2.14b) 



[See Eq. (fi~9l) .] 

Observe that S AB i s nothing but the metric tensor 
9ab- 

For symmetry class A and for given s G N, all the 
independent tensors T IJK of rank (0, 2s) of the complex 
Grassmannian can be expressed by taking the products 
of the building blocks, 

S(Aa)(bB)*{Cc){dD)»... := Z ab Y BC Z cd Y DE ' ' ' Y ZA 

(2.15a) 

made of / pairs of Z and Y, according to the decomposi- 
tion 



T, 



IK/ „ E^ = s - ( 2 - 15b ) 



For the cases of Sp(N)/XJ(N) (symmetry class C) and 
SO(2A r )/U(A 7 ") (symmetry class D) we need to replace 
Eq. (j2~T5aT) by 



S 



(i) 

...(bB)*(Cc)(dD)' 



Cc^Cc^Bc^Cd ' ' ' 



. (2.15c) 



Another simplification attached to the Kahler mani- 
folds for symmetry classes A, C, and D is that contraction 
between S"*' and the dzj's reduces to a trace, 

S tAa)(bBY(Cc)... d fi l L PAad f _ L . 1 vlB d ^ 3 L PCc ' ' ' 

= tr«A; 2 A+ •••), (2.16a) 
where we have introduced 



(2.16b) 



for the complex Grassmannian (symmetry class A). For 
Sp(JV)/U(JV) (symmetry class C) and for SO{2N)/\]{N) 
(symmetry class D) one must impose the additional con- 
straint 



(A+) T = (A-f = -ft A 



(2.16c) 



Any member T of the family of tensor fields con- 
structed in Eq. (|2.15[) has the remarkable property that 
it vanishes under the action of the covariant derivative 
on any one of its indices 



VpT 



KL- - — dpTi^ 



JKL— 

piTqjkl- 



PJ ± IQKL-- 



0. 



(2.17) 



Here, Y A BC are the components of the Levi-Civita con- 
nection whose non-vanishing entries are given by 

TV - 9 ap ' d c9p , b , T a ' b , c , = g a > d* c9pb , . (2.18) 

This is also true for the NLerM on the real Grassmannian 
0(p + q)/[0(p) x O(q)] which is related to the symplec- 
tic symmetry class of Anderson localization, and for the 
cases of simpler target manifolds such as O(N) /0(N — 1) 
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or £p n + m ~M m f or which the tensor fields T arc just 
products of the metric tensor. In all these cases, the 
fact that all high-gradient operators have vanishing co- 
variant derivatives greatly simplifies the computation of 
their one-loop RG flows. 

For a given number of derivatives, 2s, there are many 
degenerate operators of the type (|2.12[) with T IJK given 
by Eq. (|2.15b|) that share the same anomalous scaling di- 
mensions 2s at the metallic fixed point. They are spec- 
ified by the choice of } in Eq. (|2.15b[) as well as that 
of the greek indices in the d^z^s. This degeneracy is, 
however, lifted in perturbation theory as we depart from 
the metallic fixed point. 

We close this section by noting that it is more conve- 
nient to introduce the conformal indices + and — 



<9± = 9 x ±id y . 



(2.19) 



instead of the Euclidean indices [i = x,y in 2dr 



III. THE CO VARIANT BACKGROUND FIELD 
METHOD 

A. Kahler normal coordinates for NLctM 

To perform the RG program, we divide the fields (z* , z) 
in the NLctM into slow (ip*,ip), and fast (ir*,ir), modes. 
To this end, we write 



(3.1) 



where we require that (tp* , ip) satisfy the classical equa- 
tions of motion. The fast modes (tt*,tt) represent the 
quantum fluctuations around the classical background 
(ip*,ip). They are integrated out, thereby generating an 
effective action for the slow modes (ip*,ip). In practice, 
this effective action is computed order by order within an 
expansion, the so-called loop expansion, which is nothing 
but the cumulant expansion. One drawback with the sep- 
aration (|3.ip into fast and slow modes is that the pcrtur- 
bative expansion of the action in powers of the fast modes 
(it* , 7r) with respect to the action evaluated at (ip*,ip) 
need not transform covariantly under reparametrization 
of the manifold. This danger can be avoided by choosing 
Kahler normal coordinates (KNC) in terms of which we 
denote the fast modes by (ui*,uj). The KNC coordinates 
are defined by demanding the manifest covariance of the 
Kahler potential when expanded in terms of (ui* , oj) i 65 i 66 
The perturbative expansion of the action in powers of 
the KNC about the background field is then manifestly 
covariant, i.e., only covariant objects such as the met- 
ric, the Riemann tensor, the covariant derivative on the 
manifold, and so on appear in the expansion, as we shall 
shortly see. One advantage of this RG scheme, the co- 
variant background field method, is that no redundant 
operator is induced upon rcnormalization. This is a very 
useful property when dealing with the mixing of the com- 
posite operators induced by an RG transformation^ A 



related advantage is that there is no need to break the 
symmetry of the NLcrM to tame IR divergences. 

The relation between the generic parametrization 
(it*, it) and the KNC parametrization (io* ,co) of the fast 
modes is non-linear. For examples, 



i(f|/) = 9/ + fl/ 



(3.2a) 



&,,(V>*+7T*,V + 7r)=ff i .,(V>*,V) 



(3.2b) 



R t * ]kl ry li ™ ki + ■■■ 



As promised, the right-hand side is covariant under 
reparametrization as the pair (-D* , D^) denotes the holo- 
morphic and antiholomorphic covariant derivatives 



(3.3a) 



= d^+r^ b ,w* b d ti r a , (3.3b) 

respectively, with the components 

R a bcd' = ~ ^d*r a & c j R a b*c*d = ~dd,F a b* c * i (3-4) 



R ED A B = 9ec^dab an d so 011 °f the- curvature tensor. 
The expansion of the action in powers of the KNC is 



S[ip* +ir*,<ip + ir]-S[i>*,i>) = 



(p,il>)(w«w%Vd IA P 



(3.5) 



B. Canonical form of the metric tensor and 
vielbeins 

The covariant Gaussian expansion (|3.5p about an ex- 
trema (ip* , ip) of the action can still be simplified fur- 
ther by performing a local GL(M, C) transformation of 
the tangent space at (ip* , ip) that brings the metric ten- 
sor at (ip*,ip) to the canonical form. To this end, in- 
troduce the fast degree of freedom (£* , () by the linear 
transformatio n 67 ' 68 

r--=w* a a ,e a \ w* a = C\*e a \ (3.6a) 
C:=w a a e\ w « = (\e a , (3.6b) 



by demanding that 



(3.6c) 
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where the transformed metric is given by 

= diag(l Mi , ) , M 1 + M 2 = M. 



(3.6d) 



The two background dependent matrices ( Q *e a ) and 
( a ,e a ) from GL(M, C) that implement this transforma- 
tion in the antiholomorphic sector are inverse to each 
other. The same is true of the matrices ( a e a ) and ( a e a ) 
in the holomorphic sector. These four matrices are called 
the vielbeins and must obey 



9ab* 



6* 



%h-* b 'eb* & e a =9ab*, (3-7) 

by condition (|3.6cp . From now on, latin letters with a 
hat refer to the coordinates of the Kahler manifold in 
the vielbein basis (|3.6[) . Under transformation (|3.6|) the 
covariant expansion (|3.5j) becomes 

S[il>*+ir*,il> + ir]-S[il>*,il>] = 
1 



2irt 

1 
2?rf 



(3.8) 



where yet another pair of covariant derivatives 



D;C a = d^-C K b M^\ (3.9a) 

D,C = d^C-C^iA^f, (3.9b) 

have been introduced. The field is the U(M) gauge 
field defined through the U(M) spin connection a w b c . , 



by 



& e a {d c . a e b ), 
M{d ca e b )- 



(3.10a) 
(3.10b) 

(3.10c) 

[If the metric 77 is pseudo Riemannian, U(M) should be 
replaced by U(M 1 ,M 2 ).] The emergence of an U(M) 
gauge structure originates from the fact that the viel- 
beins are not uniquely specified by condition (|3.6c[) . As 
the physical result cannot depend on the specific choice 
of gauge, i.e., the specific choice for the vielbeins, the 
dependence on the spin connection should occur through 
the U(M) gauge invariant field strength 

~(F ) b = 

- & e a R r adc ,d„r*d^ d r e b + & e a R r adc ,d^ c * d^ d r e b . 

(3.10d) 



C. One-loop beta functions 

We are in position to integrate out the fast degrees of 
freedom for the complex Grassmannian, SO(2A)/U(A), 
and Sp(2AT)/U(AT). The partition function for slow and 
fast modes is 



/ DlCde- 51 **'*^'^. (3.11a) 



Z 



The action S[ip*, ip, (*, C] (= S[tp* + 7r*,V + 7r]) is sepa- 
rated into three contributions 

s[p,4>, C*, C] := S[p,4>] + SolCQ + Sj [V*, 1>, C,C], 

(3.11b) 

where S [(* , £] is the free Gaussian part for the fast 
modes, and S l [ip* , ip, £*, C] represents the interactions be- 
tween the slow and fast modes. Integration over the fast 
mode (£* , Q is performed with the help of the cumulant 
expansion 



(3.12a) 



(3.12b) 



where 



«"•»» - U^t-^' ■ (3 ' 12c) 

To evaluate the cumulant expansion, we introduce the 
Green function 



d d k e lk i x ~y) 



V & »* 27rtG (x - y), 



(3.13) 



where A = a -1 is the UV cutoff. The leading con- 
tribution to the increment of the action of the slow 
mode (ip* , ip) resulting from integrating out the fast mode 
(C* 1 C) within the second-order cumulant expansion is 

5S[r,i>] = -G o (0) J r R a , b d^»d^* a , (3-14) 



with 



R 



a * 6 



DC _ DC 

' a*c*b ' c*a*b' 



(3.15) 



the Ricci tensor. 

For symmetry classes C and D, the Ricci tensor is pro- 
portional to the metric (the target space is an Einstein 
manifold) and given by 



(3.16) 
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Thus, from Eq. (|3.14[) , the one- loop beta functions in 
(2 + e)d for the NLerM coupling (3(t) = dt/dl is 



(3{t) = -et+{M-d)t 2 +0(t 3 ). 



(3.17) 



Here, dl is the infinitesimal rcscaling of the ultraviolet 



cutoff, a ->■ ae +dl for which G o (0) = dl/2n. This one- Here, 
loop result agrees with the known three-, four-, and five- 
loop results . 35 ! 36 ' 37 Another check comes from the group 
isomorphism s] 36 ' 63 



U(2)/[U(1) x U(l)] 

U(4)/[U(1) x U(3)] 
Sp(N)/V(N) 



SO(4)/U(2) 

Sp(l)/U(l), (3.18a) 

SO(6)/U(3), (3.18b) 

SO{-2N)/V(-N). (3.18c) 



One verifies that the beta functions match if the 
couplings are appropriately rescalcd for isomor- 
phisms (|3. 18a|) and (|3.18b|) and with the additional 
substitution t — > —t for isomorphisms (|3.18c|l . The 
nontrivial zeros of the beta function in symmetry classes 
A, C, and D are 



< t* = 



0{e 2 ). 



(3.19) 



transformation law 

o\ s \z*,z) 

0\ s \z*,z) 



tdlA\fG<f\i)* : >il)). (4.1a) 



denotes all the terms that are of 



order k in ( A in the KNC expansion of 0\ 8 ' (z*,z), and 
{(••■)) = <e- Sl (--)>o/( e " Sl )o ( 4 ' lb ) 



is evaluated up to first order in tdl [recall Eq. (|3.12c[) ] . 
Once the mixing matrix 



(s«) 



(4.2a) 



is known, its eigenvalues yield the spectrum of 

anomalous scaling dimensions 



x ( ; s) = 2s-a { /h 



(4.2b) 



and the spectrum of anomalous scaling dimensions 



IV. ONE-LOOP RG FOR HIGH-GRADIENT 
OPERATORS 

A. Strategy 

We now turn to the computation of the domi- 
riant anomalous scaling dimension x\ associated to 
any high-gradient operator of type (|2.12p . where 

the collective index i runs over any allowed choice 
for in Eq. (|2.15b() and over the greek indices in 
d ^z 1 d v z J d p z K ■ ■ ■ . The one-loop anomalous scaling di- 

(s) (s) 

mension x\ of 0\ is extracted from the one-loop RG 



y\ s) :=d-x\ s > (4.2c) 

(s) 

associated to the high-gradient operator 0\ . Any eigen- 
operator of Eq. (|4.ip is relevant (marginal, irrelevant) 

when its eigenvalue yy' is positive (zero, negative). We 
are going to construct the eigenoperator of Eq. (|4.ip with 
the largest eigenvalue of the mixing matrix A^ s ^ up to 
one-loop accuracy close to the metallic fixed point. 

To this end, we need to expand 0\ (z*,z) in the 
KNC. By taking advantage of the fact that all the co- 
variant derivatives of the tensor Tjj (z*,z) vanish and 
that ([dz 1 ]^) = 0, we infer that 



([rw...^^/-]^) = (iTuJc) d^d^ ■ ■ ■ + T LL ^,r) ([d^] cl [d v z J ] Q1 )- ■• + ... . 

(4.3a) 



Thus, all we need are 



{[d^*%[d v z*%) = +tdl5 ll ,_ v d lx r k d„r l Rk* e i- i \ (4.3b) 
([VW^'W = -tdlS^d^ k d v r l R l3 'i*k, (4.3c) 
{[d^z 1 ] [d u zi] ) = +tdl5^_„d^ k d^ l R l k i h (4.3d) 
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on the one hand and 



(z*,z) 



= tdl ^2 R l * t* a * b*T bi ... brb «... b ,_ ia , b * +i ... b , (i/j*, ip) 

1=1 

r 

= tdl £ R\% T bi ... bi _ iab MI ... 6 . (V*,^), 



(4.3e) 



i=l 



on the other hand. Here, we are using the conformal 
indices defined in Eq. (|2.19[) . Note also that since holo- 
morphic V, and antiholomorphic V Jt covariant deriva- 
tives do not commute in general, we have two distinct 

representations for ([T IJK ...]^ 2 



RG transformation 



The RG transformation law (|4.3[) obeyed by any high-gradient operator for a NLcrM on the Hermitian symmetric 
spaces U(p + g)/[U(p) x \J(q)}, SO(2A)/U(A), and Sp(A)/U(A) becomes explicit with the help of the 
(a) intra-trace formulae 



tr (K] cl M[Al] cl M)) = - 2 —^tdl5^_ u [dhr (A^Af A r v M T ) + tr (A^Af) tr (AIM) 



tr([A^] (1 M[A^] cl Af 



2- tf 2 

~2-^_ v 



+tr (AIM) tr (^JSO + tf 2 tr (AJ l M T A^Af)] , 

tr (A+ r A; T AA) tr (M) - 0tr (A+ r A,; r AM/( T ) 
tftr (A+M^A^AO +tr (A/)tr (AM"^)] , 



(b) inter-trace formulae 



tdlS^_ v [tftr (Ayf* A^M) - tr (A^Af A T V M 



and (c) diagonal contributions 



2 — i9 

tdl6„ „ 

2 



tr (^A/^Af ) + tftr (^A/^Al)] , 

tr (A+ T A- T NM) - tftr ( A~ T A^ T MA/ 1 ) 
tftr (A+ T A~ T A^M) + tr (A~ T A+ T MAf)} , 



( [r 0l ... a#6 .... t . (z*,z)} \ = -tdis (Af - d) T ai ... agK .„ b . (p,1>), 



(4.4a) 
(4.4b) 



(4.5a) 
(4.5b) 

(4.6) 



where M and Af represent arbitrary rectangular matrices of the proper size and the matrix-valued A± were defined 
in Eq. (|2.16[) . This result, that follows from substituting the components of the curvature tensor (|3.4[) when expressed 
in the stereographic coordinates Q2.8[) into the master formula (|4. 3[) , agrees with the one-loop RG result for symmetry 
class A (1? = 0) in Rcf. [S ^ 



C. Diagonalization of the RG equation within the 
subspace of maximal number of switches 

As is shown in Appendix [B] the RG transforma- 
tion law (|4.3[) has an upper (or lower) triangular struc- 
ture when appropriate quantum numbers ("number of 
switches") are defined. Since the number of switches 
never increases under an RG transformation, we can 



solve the eigenvalue and eigenvector problem (|4.3p for 
each subspace with a fixed number of switches. We will 
limit ourselves to operators with the maximal number 
of switches, as the eigenscaling operator with the most 
dominant RG anomalous scaling dimension is known to 
be among them for symmetry class In this subspace, 
any operator can be expressed as a polynomial of the ob- 



13 



jects 



{tr [(A+AZ) 1 ] , tr [(A±A+) 1 ] } 



1=1,2,. 



(4.7) 



i.e., conformal indices in a single trace are completely al- 
ternating. It can also be shown that the spaces spanned 
by {tr [(A+AZ) 1 ] }, , 2 . and {tr [(AtAz) 1 ] } l=12 _ are 
not mixed by the RG transformation. We can thus con- 
sider these two subspaces separately as long as we do not 
require composite operators to be Hermitian. 

We now compute the action of the RG transforma- 
tion (|4.3[) on a composite operator of the form 



(9 



(x) 

M 



i=i 



(4.8) 



where all Qi on the right-hand side are a short-hand no- 
tation for objects of type tr [(A^Al)'] without loss of 
generality. Since the diagonal contributions (c) to the 
RG eigenvalues are trivial, we focus on the intra-trace 
contributions (a) and inter-trace contributions (b) for 
the moment. If we define the off-diagonal mixing ma- 

trrx£« = (£W {r , } ) by 



(a),(b) 



MR$ }W} 0$ } {tl,*,H (4-9) 



i.e., we have only collected the off-diagonal contributions 
(|4.4p and (|4. 5[) to the RG transformation law obeyed by 

OjA (z*, z), the action of R^ can be cast into the form 



d 



on 



2lnft 



k k 
d 2 



d 



+ (2 - tf 2 )(z + n)n t n r 



on 



l+n 



(4.10) 



The eigenvalue problem (|4.10[) can be solved by brute 
force for small s in which cases we have verified that the 
sign of the eigenvalues of the off-diagonal mixing ma- 
trix for SO(2N)/U(N) arc opposite to those for 
Sp(N)/XJ(N). In particular, the largest eigenvalue of 
i?( s ' for SO(27V)/U(iV) is given by minus the smallest 
eigenvalue for Sp(iV)/U(jV) and vice et versa. This fact 
is consistent with the group isomorphism (|3.18b[) . and 
the fact that the beta functions for the corresponding 
NLerM match upon the replacement t — > — i. 

More generally, one verifies that 



= s 
A (s) = 2s 2 



2(s 2 

2 



2s 



+ s{p + q) 

SN, (0 = 

sN, (0 ■■ 



(i? = 0), (4.11a) 
hi), (4.11b) 
-1), (4.11c) 



is an eigenvalue of the off-diagonal mixing operator R^ 
for any given s with the eigenvector 

>3 V 



M 3 



Eniw. (^ +i), (4.12b) 



M j 



2-' 



En^w. (,?=_l) - (4i2c) 



{r} 3 



As it is known that \( s ' is the largest eigenvalue of R^ 
in symmetry class A, we believe this to be also true for 
symmetry classes C and D. If so, the largest RG eigen- 
values for a given number s of gradients is given by, after 
adding the diagonal contributions (c), 

a (s) = A w -s(p + q) = 2(s 2 -s), (0 = 0), 

(4.13a) 

a w = A w - s(JV - 1) = s 2 - s, (0 = +1), 

(4.13b) 

a w = A w - «(JV + 1) = 2(s 2 - s), (<? = -l). 

(4.13c) 

Note that these results are valid up to one loop before 
taking the replica limit. Equation (|1.3[) then follows by 
combining Eq. (|3~1"5|) with Eq. (|4~2b]) . 
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APPENDIX A: GROUP-THEORY VERSUS 
GEOMETRIC REPRESENTATIONS OF THE 
NLctM 

The purpose of this appendix is to bridge the geomet- 
rical representation of the NLctM used in this paper and 
the one in terms of the so-called "Q-matrix" . 

First, recall that an element Q in U(p+g)/[U(p)xU(g)] 
can be written as^ 

Q = U- 1 AU, A = diag(I p , — 1 9 ), U 6 U(p + g).(Al) 

As noticed by Lerner and Wegner^ high-gradient oper- 
ators in the Q-matrix representation as well as the action 
of the NLerM can be conveniently expressed in terms of 
the objects 



1 



P +T (d^Q) P- T 



(A2a) 
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where t — ± and the projectors P^ 



P±:=±(l±Q), (P ± ) 2 =P ± (A2b) 

were introduced. A generic high-gradient operator is then 
built by taking a product of 

tt(A+A~A+-.-). (A3) 

For example, the action of the NLerM is 



S = 



16nt 
1 

4?rt 



* [(d.Q) (^Q)] 
tr (A; A; 



(A4) 



In order to establish explicitly the connection be- 
tween the Q-matrix representation (|A2[) - (|A4[) and the 
geometrical representation in the stereographic coordi- 
nates (|2.8|) . we parametrize P + and P~ = P + — Q in 
Eq. (|A2b|) through a, (p + q) x q matrix <& ai according 



ti 



,()!) 



P + 



af3 



E*~ 



4 = 1 



m $}i (A5a) 



where $ satisfies the constraint 

$t$ = I 9 . (A5b) 

The matrix <I> is thus a collection of q orthonormal vectors 

{$ a G C p+9 } a= i q . We can use this constraint to fix 

the upper q x q block of <&. We write 

$ = f I (A6a) 

where $ ai , <p Ijt X ij} and ^ zj X jk , are (p + q) x q, p x q, 
q x q, and p x q, matrices, respectively, and the indices 
i, I, and a run over the sets i = l,...,q, I = l,...,p, 
a = l,...,p + q, respectively. Next, we express X in 
terms of ip by making use of the constraint (|A5bj) , 

X* m Z-^X mk = 5 lk (X^Z^X = l q ). (A6b) 

Here the Hcrmitian q x q matrix Z was introduced in 
Eqs. (|2~9b)) and i[2"l?c|) . Constraint (|A^b|) is satisfied if X 
is chosen to be 



(A6c) 

3]) in the Q-matrix 



A f = X = Z 1/2 . 

A generic high-gradient operators 
representation is thus expressed in terms of the stere- 
ographic coordinate (ip*,ip) from Eq. (|2.8[) . For exam- 
ple, one verifies that the action of the NLerM (|A4[) re- 
duces to Eq. (|2.5p with the metric (|2.9ap derived from 
the Kahler potential (|2.8p once the parametrization (|A5|) 
with Eq. (|A6[) is used. Similarly, 



tr (A+A-A+A-) = 
tr [Z (5^t) y (d v <p) Z (d^) Y (dM ] . (A7) 

Observe that the alternating structure in A + and A~ in 
the high-gradient operators has a one-to-one correspon- 
dence to that in Z and Y, 



APPENDIX B: UPPER TRIANGULAR 
STRUCTURE OF THE RG TRANSFORMATION 

We prove in this appendix that the RG transformation 
law obeyed by high-gradient operators has an upper (or 
lower) triangular structure when appropriate "quantum" 
numbers ("number of switches") are defined. A trace 
made of an even product of As, 



tr (A+ T A~ T A+ T A 

\ £71 (72 <^3 



— r 
cr 4 



(Bl) 



is said to have n T a number of switches from A^ to A_Z 
when, dividing the sequence into pairs according to 

tr [{Atl A':) {At: A-:)-.-}, (B2) 

there are n T a pairs {A^AZ 1 ^) in the sequence. The number 
does not depend on the cyclic permutations of the 
pairs (AZ.AZ&)- For a given trace of the form (|B2|) we 
can define 



E < 

T,<7 — ± 



(B3) 



We are going to show that, for the symmetry classes 
A, C, and D, the total number 



tr 



(B4) 



of switches in a composite operator of type (|2.12p never 
increases under the RG transformation law (|4.3[) . Here, 
the summation extends over the product over all traces 
of the form (|B2|) making up the high-gradient operator. 
In other words, the RG transformation law (|4.3p is tri- 
angular with respect to the quantum number (IB4|1 . 



1. Symmetry class A 

We begin with symmetry class A. This case ($ = 0) 
was already treated in Refs. [3^[3^j43| and is needed for 
symmetry classes C (# = — 1) and D (fl = +1). Consider 
any string L-^A^^ AZj R x and any string L 2 A^^ A~7 R 2 of 
4 matrices each that enter a high-gradient operator. 

First, we assume that each string belongs to two dis- 
tinct traces. We seek the one-loop RG transformation 
law obeyed by the number of switches 



contained in 



(B5a) 



tr \L x Al x A-7RA tr \L 2 A^ A~7R 2 ) . (B5b) 
By making use of Eq. (14. 5|) with ■& = we deduce that the 
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one-loop RG transformation law obeyed by Eq. (|B5b[) is 



tr 
tdl 



-tdlS. 



2^2 



A T a A~7 R X L X tr AZ A~7 R 9 L 

6 <Tl,-<T2 + 6 <tx,-<t' 2 + S a[,-a 2 ~ 5 <j[ -c 

x tr [AtlA-rR^AtlA-jR^ 



+tdl8 n „, tr f R 9 L 9 At T A~7 A± T A~7 R, L, 
+tdlS, „ tr ( A± T A~7 A± T A~7 R 2 UR, L-, 

°i>— CT 2 V 1 "i 2 "a 



tdU^tr yA+ T A~7 R 2 L 2 A~t. T A~7 R X L X 



(B6a) 
(B6b) 
(B6c) 
(B6d) 
(B6e) 



Line (|B6a[) contains 5^ _ a , + S a 
Lines (|B6bl) to jB6c) contain o„ 



, switches. 

i 

+ S 

a 2 ,—<T 2 

switches. Hence, only line (|B6a[) has the potential 
to increase the initial number (|B5a[) of switches. 
However, an increase over the initial number (|B5ap 
of switches is not possible in view of the iden- 
tity -<5 CTij _ CT2 + S ai _ a , 2 + ^ - <S CT ,_ CT , = 
2a 1 a 2 S a _ a , S a _ a , . Line (|B6a[) only contributes if 

a[ = —a x and a' 2 = — a 2 both hold, i.e., when the 
initial number of switches (|B5a[) is maximum (2). 
Consequently, the number of switches 25 a CT induced 
by line (|B6a[) equals the initial number of switches 2 
when a x = a 2 while it decreases otherwise. 

Second, we assume that each string belongs to the 
same trace. We seek the one-loop RG transformation 
law obeyed by the number of switches 



contained in 



(B7a) 



(B7b) 



By making use of Eq. (|4.4|) with fl = one verifies again 
that the number of switches n T a cannot increase after in- 
tegrating over the fast modes up to one loop. 

We have proven that the RG flows of high-gradient op- 
erators made of any given number of switches decouple 
from the RG flows of high-gradient operators made of 
a lesser number of switches in symmetry class A. As a 
corollary, the mixing matrix (|4.ip is triangular in sym- 
metry class A. 



2. Symmetry classes C and D 

We turn our attention to symmetry classes = 1, 
i.e., SO(2AO/U(7V) and Sp(A)/U(A). We first consider 
the inter-trace formula of the RG transformation 
applied to 



where A4 :— R X L X and Af :~ R 2 L X . It is sufficient to 
consider the contributions proportional to $ in Eq. (|4.5p 
as the remaining terms can be treated along the lines 
of symmetry class A. There are four such contributions 
(a) - (d): 

K^{K^ - K^YAAtl^ A-7AXA-7M), 

(B9a) 

<^,- CT ;(<W 2 - 5 au _ a2 )tMtI^A-7At;_A-7M\ 

(B9b) 

S„ „, (6 a , , — 6 a , ,)ti-(A-7A+ T Af T A-7MA+ T ) 7 

<Tl, — <T 1 \ ^ii^^ <T 11~ <T 2^ ^ °2 2 °1 1 

(B9c) 

S„ n , (S„, n , - S„, A± r N T A~7 MA± T ). 

(B9d) 

We assign to 

A := tr (At7A-7M) x tr (A+TA^Af) (BIO) 

the quartet of numbers 

(nf, nf, n7, nZ) (Bll) 

where nZ counts the total number of switches of type 
A T a A~77a in A. In symmetry classes D and C, we can com- 
bine trO = trO T and the symmetry {A±) T = -tiA^, to 
infer that 

(n^, nf, nZ, nZ) = (nZ, nZ, nf, nf). (B12) 

As opposed to symmetry class A, individual switches nZ 
are not separately conserved but the total number of 
switches 



= 2_. n "-r 

I7.T 



(B13) 



will be shown to be conserved. 

Without loss of generality we choose r = + in A. The 
total number of switches nf and nf are 

n% = n\ [AX A~, }+n+[M} + n% [A+ A~, ] + nf[N], 

(B14a) 

nf = nf [A+ A~, }+nf[M} + nf [A+ ] + nf [Af] , 

(B14b) 

respectively. The total number of switches nZ and nZ 
follow from using A = tr (A^MA+J x tt(A~,MA+ 3 ) 
and are given by 



n + = n + [A, MA^ ] + n + \A,M A+ ] 



n~ [A- M A+ 1 + n- [A- AfA± 1 , 



(B14c) 



(B14d) 



• • • X tr [AT. A~7M x tr ( AZ A~7Af) x 



(B8) 

respectively. 
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The number of switches (|B14I) in A changes under the 
RG transformation (|B9|) as follows. The total number of 
switches nt and nZ from line (|B9aj) is given by 

n+[(a)] = n+lA^A-J+n+lA^A^+n+lM] 
= nZ [A-,MAX } + n\ [A+ A~, ] + n+ [X] , 

(B15a) 

n±[(o)] = nt[<AA T A;,]+nt[A+A;,]+nt[A4] ' 
= n- + [A-,NAX } + nt [A+ t A~, ] + n± [M] , 

(B15b) 

respectively, where we again made use of the symmetry 
condition {A^) T = —flA^, and the fact that the num- 
ber of A^ in a trace is always even. Similarly the total 
number of switches nZ and nZ from line (|B9a[) is given 
by 

n-[(a)} = nllAf^+nllA-A+A+nllA-MAX] 
= nt[N]+nl[A-,AtA+n+[A-,MAt 2 ], 

(B15c) 

?il[(a)] = (as above with o — ► — <r), (B15d) 

respectively. Thus, the total number of switches changes 
from 

n = Y, {<\<K { ] + »t[M] + nt[Ai 2 A-,) 

a 

+n+ [M] + n- [A-,MA+ } + n~ [A~,AfA+ 2 } } 

(B16) 

to 

n[(a)} = ]T {nZAA-,MAt 2 ] + n+[<^] + n+[M] 

(T 

+n+ a [Af] + n-[A-,AX] + n~[AZ.MA+ a ]} 

(B17) 

for the contribution from line (|B9a|) . The net change 
An| (a ) = n - n[(a)] is 

An| (a) = £{<[^^]- n ;[^<] 

(7 

(B18) 

Whcn ^2,-4^1,^ = h wc find 

(T 

+n-[A-,MAt 1 ]-n-[A-,MAt 1 ]} 
= (B19) 



with the help of the prcfactor 6 a2 _^ (S au(T3 - _<r 2 )- 
Whcn ^2,-4^!,-^ = wc find 

An| (o) = £{n+[A+A;,]-n"L4;,A+] 

(7 

> 0. (B20) 

We conclude that the total number of switches docs not 
increase under the RG transformation (|B9a|) . 

Repeating the same analysis for the contributions (b) — 
(d) in Eq. (|B9I) yields that the change in the total number 
of switches is > 0. 

We leave it to the reader to verify that the same con- 
clusion holds for the intra-trace formula ()4.4|) . 

This completes the proof for the upper triangular 
structure of the RG equation with respect to the total 
number of switches in symmetry classes D and C. 



APPENDIX C: RG FLOWS OF 
HIGH-GRADIENT OPERATORS FOR NLcrM ON 
SUPERMANIFOLDS 

We derive in this appendix two "master formulae" 
needed to perform the one-loop RG analysis of high- 
gradient operators for NLcrM defined on an arbitrary Rie- 
mannian or Kahlerian supcrmanifold. Master formulae 
(|C19|) and (|C40[) are supersymmetric versions of Eq. (|4.3p 
and can be used for any target supermanifold once the 
geometrical data of the manifold are known. An expan- 
sion in terms of normal coordinate of tensor fields plays 
again an essential role. For the non-supcrsymmctric Rie- 
mannian case, an expansion in terms of Ricmann normal 
coordinate (RNC) based on the notion of geodesies was 
used for the computation of the beta function in all gener- 
alit y 32 i 57 and the RG analysis of high-gradient operators 
for 'the special case of the 0(N)/0(N - 1) NLaMii^. 
The RNC expansion can also be applied to Hermitian 
manifolds. For the (non-supersymmetric) Kahlerian case, 
however, one can also use the Kahler potential to find 
a normal coordinate system, the Kahler normal coordi- 
nates (KNC) introduced in Refs. [65l and [6(1 In this Ap- 
pendix, we will generalize these normal coordinate expan- 
sions to the case of supersymmetric target manifolds. All 
the results that follow reduce to the non-supersymmetric 
ones if we ignore the commuting-anticommuting nature 
of the objects and fermionic sign factors. 



1. NLctM on Riemannian supermanifolds 

a. Notations and conventions for Riemannian super- 
manifolds We begin with a quick review of concepts and 
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notation. We refer the reader to Ref. \5£ 
on the geometry of supcrmanifolds. 



for a textbook The supcrsymmctry of the metric tensor field implies 



A tensor field of rank (r, s) on a supermanifold has 
r contravariant (upper) and s covariant (lower) in- 
dices, T ab c d . On a supermanifold, components 
of tensor fields can be either commuting or anti- 
commuting supcrnumbers. 

In order to define their mutual statistics, the grade 
e(0) = 0, 1 is assigned to each object on a su- 
permanifold, where can be a supernumber, a 
supcrtensor field, a component of a supcrtensor 
field, etc. When the ordering of two objects Oi 
and 02 is exchanged, a factor (— l) e (°i)e(02) arises, 
0102 = (-l) e(c,l)e(O2) 2 0i. In order to define the 
grade for the component of a tensor field of arbi- 
trary rank, we first assign a grade to each of its in- 
dex. If index a corresponds to a commuting compo- 
nent of the coordinates of a point on the supermani- 
fold e(a) = 0, while e(a) = 1 otherwise. For simplic- 
ity, a statistical factor (— l) e ( a ) e ( b ) is often abbrevi- 
ated by (-l) Qb , e.g., X a Y b = {-l) ab Y b X a . The 
grade of the component T afc c d ... of the tensor field 
T is then e(T ab c d ...) = e(o)+e(6)+e(c)+e(d)+- • • & 

It is convenient to introduce a shifting rule for the 
left most index of tensors, 



arpb 

Tb 



rpab... 

■— 1 C...1 

:= (-l)°T a 6 



(Cla) 
(Clb) 



• As with all other objects on a supermanifold, 
derivatives also carry a grade. Thus, derivatives 

from the left d a and the right d a need to be dis- 
tinguished. They are related by 



~daf = (-l) ae(/)+ °/K, 



(C2) 



where / is a function over the manifold. In this 
appendix we will mainly use derivatives d a and 
covariant derivatives V a that act from the right. 



• Super transposition of matrices, denoted by 
defined by 

3tf~3 = (_1)«<JVJ. 



IS 



(C3a) 
(C3b) 
(C3c) 



For illustration, we shall apply these rules to the com- 
ponents g ab of the metric tensor field, and components 
9 

shifting rule 



b of the inverse of the metric tensor field. From the 



V = 9 a \ 
«9 c c 9 b = j\ 



a9 b = (-l) a 9 ab , (C4a) 



9 c9b 



a s h . 



(C4b) 



a 9b = (- l ) a 9ab 

= (-l) a+ab 9 b 



b9a-> 



= (-l) ab g ba 
= (-l) abb g a . 



(C5a) 



(C5b) 



An index of a component X a of a supervector field can 
be raised, lowered, and shifted according to the rules 



a X = a9b 

a X = a g b b X, 



a X = (-irX a , (C6a) 



X a — X b g a , 
X a = X b b g a . 



(C6b) 
(C6c) 



b. The NLaM on a Riemannian supermanifold The 
NLctM on a Riemannian supermanifold is defined by the 
partition function 



Z := / V[<j)}e 



1 

4?rt 



dA a a 9 b {4>) d**, 



(C7a) 
(C7b) 



where (<j) a {r)) € R* 1 x Kf represents the M + N compo- 
nents of the coordinates of point r on the target manifold, 
with A I refcring to the number of commuting coordinates 
whereas N that of the number of anti-commuting ones,— 
and t is the NLaM coupling constant. 

c. The covariant background field method and the su- 
per Riemann normal coordinate expansion In the back- 
ground field method, the covariant vector field <p a (r) is 
separated into two parts, 



(C8) 



whereby ip a (r) is assumed to be a slowly varying solution 
to the classical equations of motion which transforms like 
a contravariant vector and 7r a (r) represents fluctuations 
around the slow degrees of freedom i[i a (r). 

Two steps are needed to perform the one-loop RG pro- 
gram in a covariant fashion. First, we need to trade the 
expansion in terms of 7r°, which is not covariant, for that 
of fields £ a that transform like contravariant vectors. Sec- 
ond, we need to rotate the RNC by a viclbcin to evaluate 
the functional integral over the fast modes. 

The RNC is based on the notion of geodesies on Rie- 
mannian manifolds. Since almost all the notions on a 
Riemannian manifold such as forms, connection, the Rie- 
mann tensor, geodesies, etc., have their counterpart on 
a Riemannian supermanifold^, it is straightforward to 
develop the super extension of the RNC expansion. The 
expansion of a covariant tensor field T in terms of the 
super RNC £° is given by 
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T K ... h {i> + tt) = T v .. 6i (V) + r 6i .„ ti V c (^) r + 
1 ' 



1 



2! 

+b)(b i+1 +—+ 



6, Vcj V C2 



i=l 



')(i? b biCiC2 + (-l) b ' Cl i? b CibiC2 ) 



+ 



The RNC expansion for 9 M (?/> a + 7r a ) is given by 



(C9) 



(CIO) 



where Z)„ is the right covariant derivative for £ a 



c c b 



(Cll) 



and r a bc is a component of the connection. For a Rieman- 
nian supcrmanifold the connection can be derived from 
the metric. With the covariant expansion Eqs. (|C9[) and 
(|C10p the action is expanded in terms of C according to 



S[V,7r]-5[V] 

r \e a9b w b (D^ L + R abcd d^ d c e d^ a 



Ant 



(C12) 



d. Vielbeins on a Riemannian supermanifold To in- 
tegrate over the fast modes in the expansion (|C12|) . we 
rotate C to C with the linear transformation 





r = C Q a e a W, (C13a) 


d C:= a e a (^) a C, 


^ = b e g (V') S C, (Cl3b) 


where 




e h ~e h = S b 


a e a a e b = a <5 b , (C14a) 




a e a a e 6 = % (C14b) 


(_l)4C4+«) a g« = «g a> 


(-l) a (»+ a ) a e a = a e a . 



(C14c) 

Here, the vielbeins (and their inverse) e are defined by 



&e a a 9b b h= a e a a% b eb = a9b, (C15a) 



with 



£> M 



D 



N 



(C15b) 



D M := diag -1, • • • , -1, +1, ■ • • , +1 ,(C15c) 



D N := 



(-° i y)--(- i^)>'») 



r 



Observe that we are allowing the bosonic part of rj to be 
pseudo Riemannian, i.e., r\ is allowed to have —1 as well 
as +1 as its diagonal elements. From now on, latin letters 
with a hat refer to the coordinates of the Riemannian 
manifold in the vielbein basis. 

The covariant expansion of the action in terms of C a is 
given by 

S[iP,n]-S[ij] = 

+ • • • , (C16a) 
where yet another covariant derivative from the right 

:= 9/C + a 0yV0) S S C, (C16b) 



(C16c) 



has been introduced together with the OSp(M|iV) spin 
connection [if the bosonic part of 77 is pseudo Riemannian, 
OSp(M|iV) should be replaced by OSp(Mi, M 2 \N) where 
M1/2 is the number of ±1 in the diagonal part of 77 with 
Mi + M 2 = M] , 



(-iy» ( a i- h d c ) +r a cb % 



(C16d) 



the OSp(M|7V) gauge field 



N/2 



& (A,) h := (-lf-^d^, (C16e) 
& (A,y b = -(-lp+^^f, (C16f) 

and the OSp(M|A r ) field strength tensor 

+ & {A ll ) £ & {A v y b -\A v ) £ & {A^\ 
= -%R\cdd^ d W%- (C16g) 

The components of the Riemann tensor are here R a bcd - 
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e. One-loop beta function The integration over the 
fast modes £ is performed with the help of the cumulant 
expansion. The effective action for the slow modes ip is, 
up to one-loop, 

Z:= Jv[i;]e- S ^ +5S ^\ (C17a) 

ss[i>] = \g (o) Jdrf (-1)% d„ V,(ci7b) 

where 

Ra ■■= {-i) k(l+1)+l g kl Ruk 3 (Cl7c) 



is the Ricci tensor. 

/. RG flows of high- gradient operators The covari- 
ant expansion in terms of the super RNC applies to the 
high-gradient operators 

Tijk.. .{<!>) d^ftd^ d p <j) k ■ ■ ■ (C18) 

as well and can be used to compute their anomalous scal- 
ing dimensions. We shall assume that all the covariant 
derivatives of the tensor field T i j k (<fi) vanish. If so wc 
infer the RG transformation law 



Thus, all we need are 



[d v <tt] c ) ••• + ••• • 



on the one hand and 



T, 



(0) 



2! 3 2^ 1 b 1 -b i _ 1 bb i+1 ...b l { J-J 

I 



R" 



(-lf^R b c , c , 



(C19a) 

(C19b) 
(C19c) 

^9°\ (C19d) 



on the other hand. Here, we are using the conformal 
indices defined in Eq. (|2.19p . 



2. NLctM on Kahler supermanifolds 

g. Notations and conventions for Hermitian super- 
manifolds Some additional notation is needed to deal 
with Hermitian and Kahlerian manifolds. 

• With the convention that complex conjugation 
of a composite objects reverts the ordering, 



plex conjugation of a tensor field is given by 



T 



A x ... 



(C20b) 



v?ith5 r+s (A) = Elu=i,t<u A t A v 
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• The grade for an antiholomorphic index is equal to 
that of its holomorphic counterpart, e(a*) = e(a). 



For illustration, we shall apply these rules to the com- 
onents a g b „ of the metric tensor field. Combining the 
rule (|C20b[) with the supersymmetry of the metric (|C5[) . 

(-ir +b+ab ( a ,g b )\ a ,g b = 



a9b* 



b9a* 



(010203 • • •)* - ' ' ' (03)*(02)*(0i)*, it is useful ppne nts a g b » of the metric tensor field. Combining the 
to introduce sign factors associated to the opera- 
tion of complex conjugation. For covariant and con- 
travariant holomorphic and anti-holomorphic vec- 
tors, the action of the complex conjugation is de- 
fined by ^ ' 

(v a y = v* a , (v a y = {-iyv*. (C20a) 

For tensor fields of interest in this paper, the com- 



we infer that the metric tensor field is superhermitian 

9b>a = (9a*b)*. (C22) 
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h. The NLaM on a Hermitian supermanifold The 
NLcM on a Hermitian supermanifold is defined by the 
partition function 



Z 

S[z*,z] 



1 

2nt 



V[z*,z]e- S[z%z \ (C23a) 
d fi z*\.g b (z*,z)d fi b z, 



(C23b) 



where {z* a , z a ) represents coordinates on the complex su- 
permanifold. For a Kahlerian supermanifold, the metric 
is derived from the Kahler potential K (z* , z), 



9ab* = K(z*,z)d a d b . 



(C24) 



[Here, observe that the shifting rule for indices (|Clb[) and 
the relation between the left- and right-derivatives (|C2[) 



are compatible.] 

i. The covariant background field method and the 
super Kahler normal coordinate expansion In the 
background field method, the covariant vector field 
(z* a , z a ){r) is separated into two parts, 

z* a = i>* a + 7r* a , z a = V a +7r a , (C25) 

whereby (ip* a , ip a )(r) is assumed to be a slowly varying 
solution to the classical equations of motion which trans- 
forms like a contravariant vector and (7r* a , 7r a )(r) rep- 
resents fluctuations around the slow degrees of freedom 

The expansion in terms of (it* , it) is traded off for one in 
terms of normal coordinates (w* , w), the super version of 
the KNC. The super KNC expansion for a tensor field of 
rank (0, r + s) with r holomorphic and s antiholomorphic 
indices T bi ... brb ,... b ,(z* , z) is 



bi ... Kcl ... c ,(z* 7 z) = T bi ... Kc .... ci (r,^)+T bi ..^^ 

^...6_c?... c! V fel V foW fe2 W fel + l -T bi ... Kcl ... c , V,. V liW * h w* h +T b ^.. bc ,... c ,V kl V l; w 



r 

V^(-l)(°+ b *)( 6 *+i+'"+ c »)T i 



b 1 ---b i _ 1 ab i + 1 ■■■b r c 



1—°% h i k i l l 



w* h w kl 



T, 



, i ... MI ... c .(^,^)+r bi ... brCr ... c; V fci «; fcl +T h ... brC ,... c ,y i . i 



w" 



+ \T bi ... Kcl ... c *y kl V ko w k >w^ + \T b ^ +T ( 



2 b 1 ---b r c*---c* 



l bi ... brC *... c * V|;V fel 



» V; f Vi,w 



fcl *fl 



,...^i? a * c . irfc w 



fcl„„*il 



E(— 1 U°+ c 0( c i+H I" c *)t 
I ^ - t & 1 ---& r c*---c*_ 1 ( 

i=l 

The super KNC expansion for d ll ('tp A + w A ) is given by 



(C26) 



(C27) 



where D is the right covariant derivative for (w* j ,w 3 ) is expanded in terms of (w* a , w a ) according 



to 



= d^vP +T\ b w b d^ a , (C28a) 



*w*D* u = iw*V a «d^* a 



= 9>* +T^ a , b ,w* b d^* a 



(C28b) 



with T A BC the coefficient of the connection that satisfies 

(r a 



S[lp* + 7T* , ijj + 7r] — S"^* , V] = 

1 



(C29) 



2^c+6+a(6+c)+bc r a *f,» c » [see (|C20b|) ]. With j. Vielbeins on a Hermitian supermanifold Intro- 

the covariant expansion Eqs. (|C26[) and (|C27[) the action duce the fast degrees of freedom (£*,C) by the linear 
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transformations 





C -=w a a e a , 


u. a = C a a e a , (C30a) 




& C-= & e a a w, 


b u; = b e, b C, (C30b) 


and 








C* a := w* a a ,e a *, 






:= a *e „. a w*, 


b u>* = b "e t b (* (C30d) 


where 








P& pb _ sb 
a° a° a u i 


a e a e — a i \\>OUe) 






Ss Oi _ ax fT"}nn 
a b ~ b' l,^<5UlJ 


(- 


_]\a(a+a)_ go _ a-_ 


(-l)«( a + a ) a e a = a e a , 






(C30g) 




= (-l)( a+a )° a .e a *, 


( a e a )* = (-l)( a+a >Ve a *. 






(C30h) 



These definitions of the viclbcins are equivalent to the 
non-supersymmetic cases except the last two lines. The 
vielbeins (and their inverse) e are defined such that 

w\g b . b w* = C & V- b . h C (C30i) 

where the canonical form of the metric r\ is given by 

a%, = diag(l Ml , -I M2 ,I N ). (C30j) 

Note that we are allowing the bosonic part of r] to be 
pseudo Riemannian, i.e., 77 is allowed to have —1 as well 
as +1 as its diagonal elements. From now on, latin let- 
ters with a hat refer to the coordinates of the Hermitian 
manifold in the viclbcin basis. 

The covariant expansion of the action in terms of 
(C\ C a ) is given by 

S[ip* +ir*,ip + ir}-S[<ip*,i>) = 
1 

2~7rf 

c*V<)(cv*) 



Ri* 



i'*jkl* 



2lTt 

J(C i "-,e i *)9 M ^V J '(rS»e i *) 



2 
1 

where yet another covariant derivative from the right 



(C32a) 
(C32b) 
(C32c) 

^C* a -C*V(^) S * (C32d) 
has been introduced together with the U(Mi,M2|iV) 
gauge field A,. The gauge field with different index- 
structures are related by 

lp+o^f, 



D*X* a 



a (D, = CD, =► a (A,) b = - 



(C33a) 
(C33b) 



The gauge field is skew superhermitian, 

~ b {A-n)c' = — (c(A") *) 



(C34) 



where indices of A, are raised and lowered by r\. The 
explicit form of the gauge field in terms of the vielbeins 
is given by the U(Mi, M%\N) spin connection 



a 6*c* 

whereby 



a e Q ( a e b d c ,), (C35a) 

a e a [( a e b 5 c ) - (_i)»+»H-(H-S)c T \ c r e b ] , 

(C35b) 

(-l) c ( b + a V a . c )*, (C35c) 



& {A»t =■ au b c *d,r c + a" b c d^ c - (C36) 
The U(Mi, M2|iV) field strength tensor is given by 

a(^) S -^ a (^) b -9, a (^) b 

+a(Afj.) C c(Au) b ~ a{A v ) C c{A,) b 

= -(-l) a+ar h e a R\ dc , 8„r c d^ d r e b 

+ { _ ir+ ar a g„ ^.c ^ ^ ^ 

k. One-loop beta function The integration over the 
fast modes (£* , C) is performed with the help of the cumu- 
lant expansion. The effective action for the slow modes 
(ip*,ip) is, up to one-loop, 



(C31) 

where 



Z := J T>W,il>] e S[r A>]+ss[r M ^ ( C38a ) 
5S[r,^] = -G (0) [ R a , b d^ b d^* a ,(CZ8b) 



(C38c) 



is the Ricci tensor. 

I. RG flows of high-gradient operators The covari- 
ant expansion in terms of the super KNC applies to the 
high-gradient operators 



T IJK ,Sz\z)dz I tl dz J v dzf 



(C39) 



as well and can be used to compute their anomalous scal- 
ing dimensions. We shall assume that all the covariant 
derivatives of the tensor field T IJK (z*,z) vanish. If so 
we infer the RG transformation law 
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[T IJ Jz*,z)d ll z%z J --] ) = ([T„...(^^)] c2 )V J 5 [/ ^-..+T„...(^^)([V / ] cl [d„z J ] 



Thus, all we need are 



= -^ ltt _ v R* be . d a lt fi> d a v r cb 9 jm , 

([d„z H ] (1 [d^] (1 ) = -ttB5„ i _ v nr Vcd .d ll r%r h '9 j , 



on the one hand and 



(C40a) 

(C40b) 
(C40c) 
(C40d) 
(C40e) 



T bi ... brC .... c .(z*,z)\ ) = -*£fiX;(-l) ( ^ 6<)(64+1+ "^ )r 6 1 ...6 i _ 1 a6 <+1 ... ir c J ...c:i2 6 t * llJ * ,1 ff fcl 

s 

= -^E(- 1 ) (a+Cl)(C!+1 + ''' +Cs)T ^---^ r .. c: _ 1 a.c :+1 .. < fi Q *c*q fcl fcl 3* il , (C40f) 



on the other hand. Here, we are using the conformal from the RNC expansion as well. A subtle thing here 



indices defined in Eq. (|2.19D . 

Since a Kahlerian supermanifold is a Riemannian su- 



is that the RNC and KNC expansion for a tensor field 
Tijk...( z * \ z ) an d d t _ l z A are different However, the 



permanifold, one can make use of either the RNC or KNC RNC and KNC expansions are identical when applied to 
expansion. Indeed, all the results in this appendix de- 



the high-gradient operators 
rived from the KNC expansion can actually be obtained same results. 



and thus give us the 
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